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Abstract: We elaborate on a recently found SU{5) relation confined to the up-(s)quark 
flavour space, that remains immune to large quantum corrections up to the TeV scale. We 
investigate the possibilities opened by this new window on the GUT scale in order to find 
TeV-scale SU{5) tests realizable at the LHC. These SU{5) tests appear as relations among 
observables involving either flavour violation or chirality flip in the up-(s)quark sector. 
The power of these tests is systematically evaluated using a frequentist, p-value based 
criterion. SU{5) tests in the Heavy supersymmetry (SUSY), Natural supersymmetry and 
Top-charm supersymmetry spectra are investigated. The latter scenario features light stops 
and scharms and is well-motivated from various hve-dimensional constructions. A variety 
of SU{5) tests is obtained, involving techniques of top polarimetry, charm-tagging, or Higgs 
detection from SUSY cascade decays. We find that 0(10) to 0(100) events are needed to 
obtain 50% of relative precision at 3cr significance for all of these tests. In addition, we 
propose a set of precision measurements in ultraperipheral collisions in order to search for 
the flavour-changing dipole operators of Heavy supersymmetry. 
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1 Introduction 


The Standard Model (SM) of particle physics describes to a rather good accuracy all kinds 
of elementary matter observed up to now - apart maybe from the elusive dark matter 
present in our Universe. A fairly fascinating feature of the SM is that these matter fields 
fit into complete representations of the SU{5) gauge group [1]. More precisely, the fields 
{Qi, Ui, Ei} and {Li, Di} can be respectively embedded into three copies of the 10 and 
the 5 representations of SU{5) This feature can be interpreted as a hint that the SM 
gauge group SU (3) x SU (2) x [/(I) arises as the low-energy limit of a Grand Unified Theory 

(GUT) containing the SU{5) sub-group at some higher scale (see Refs. [2, 3] for reviews) 
2 

Unravelling whether or not Nature is 5?7(5)-symmetric at short distance constitutes 
a challenging open problem of particle physics. Many realizations of (S'?7(5)-like GUTs 
are possible, with various low-energy consequences. Moreover, the quantum corrections 
between GUT and low-energy scales potentially hide a lot of information. A SU{5) test 
has thus to face both challenges of avoiding a large model-dependence and being doable in 
the real world. 

Besides the equality of the gauge couplings, a consequence of SU (5)-like unification of 
the matter fields is the relation 

Vd = y\ (IT) 

between the Yukawa matrices of the down-type quarks and charged leptons. The fact that 
we do not observe equal masses for leptons and down-type quarks at low energy may be 
explained by the renormalization group (RG) flow, which does not respect SU{5) symmetry 
below the breaking scale. Moreover, Eq. (1.1) is also modified by finite corrections at the 
GUT scale, coming in when integrating out the heavy degrees of freedom of the theory. It 
is thus necessary to study to which extent this relation can be fulfilled at both high and 
low energy, and a quantity of work has been done on the subject [5-32]. This is a typical 
example of SU (5) test where observables are simple but model-dependence is large. 

In the hypothesis that supersymmetry (SUSY) is present at a somewhat low scale, new 
forms of matter exist, that potentially offer new insights on SU (5) symmetry. We will work 
within this hypothesis throughout this Paper. Supersymmetry is a natural companion of 
SU (5) unification as the simplest supersymmetric models (such as MSSM, NMSSM, ...) 
lead to gauge coupling unification with a good precision. This happens when no new degree 
of freedom shows up between the electroweak (EW) and the GUT scales, but also if new 
spacetime structure appears at some intermediate scale, or if additional fields in SU (5)-like 
representations are present. Apart from that, SUSY also naturally stabilizes the EW scale 
and provides candidates for cold dark matter in our Universe. 

In the SU (5) framework, the two Higgs supermultiplets of the MSSM, denoted and 
Hu, need to be embedded in a 5 and a 5 representation, denoted Hd and Hu, respectively. 

^Note that this is not the case for the Higgs doublet(s), for which a splitting mechanism is required. 

^Some early SU{5) realizations are exclnded by proton decay [4]. 
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Interactions between matter and Higgs fields are then given by the superpotential 

W = + Al^TiulOilOj . (1.2) 

After SU{5) breaking, assuming that unwanted Higgs triplets are heavy due to some mech¬ 
anism, the superpotential reads 

W = vHHdQiDj + yfHaL.Ej + y^HuQiUj . (1.3) 

At the GUT scale, Eq. (1.1) is still verified, with yu = y\ = ^u- Furthermore, supersymme¬ 
try needs to be broken, such that the Lagrangian also contains the SUSY-breaking scalar 
trilinear terms, 

>Csoft ^ auhuQU + Qdhdqd + a^hdie , (1.4) 

and SUSY-breaking scalar masses, 

Csoit D -TMi^q-u*M^u-d*Mld-e*M^e-tMli, (1.5) 

where q, u, d, i, and e are the supersymmetric partners of the SM fermions. With the 

assumption that the source of SUSY breaking is S'[/(5)-symmetric, SUSY-breaking terms 
can be parametrized by a singlet spurion and the aforementioned SU{5) relation is also 
verified for the corresponding trilinear terms, 

ad'^a\. (1.6) 

It would be again tempting to use these relations in order to test the SU (5) GUT hypothe¬ 
sis. However, contrary to SUSY-preserving parameters, it turns out that the RG evolution 
of Eq. (1.6) is much more model-dependent than the one of Eq. (1.1). It depends on the 
other parameters of the model, on the hypothesis about SUSY breaking, on the hypothesis 
of having no new degrees of freedom, and so on. It may even be dominated by the SUSY 
breaking hidden sector dynamics. Therefore this relation can hardly constitute a generic 
test of the SU{5) GUT hypothesis. The scalar masses satisfy the SU{5) relations 

(1.7) 

The same arguments from model-dependence apply to these relations as for the one of Eq. 

( 1 . 6 ). 

Testing the SU (5) hypothesis in a way as model-independent as possible, even with low- 
energy supersymmetry, seems therefore to be a rather arduous task. However it turns out 
that an additional 51/(5) relation, first identified in Ref. [33], opens intriguing possibilities. 
Building on this observation, the primary goal of this work is to find realistic tests of the 
51/(5) hypothesis, and to evaluate their potential for the Large Hadron Gollider (LHC) 
era. Let us stress that our study opens a new possibility with respect to those relying 
on the quark-lepton complementarity of Eqs. (1.6) and (1.7) (see, e.g.. Refs. [34-36]). 
Such correlations, although certainly interesting in specific models, can be hardly used as a 
generic test of the 51/(5) hypothesis, as the RG corrections received by quark and leptons 
are fundamentally different. 


- 3 - 


The paper is organized as follows. In Sec. 2, we review the sector of up-type squarks 
and discuss the SU{5) property which is to be exploited in the subsequent Sections. The 
persistence of this property is closely examined and quantified in Sec. 3. The relation is 
found to be exact up to 0(1%) relative discrepancies. Sec. 4 comprises the strategy and 
tools needed to setup SU{5) tests at low energy. It includes in particular the up-squark 
effective Lagrangian, and a systematic frequentist method to evaluate the power of the 
proposed SU{5) tests. 

We then investigate various SUSY scenarios. In Sec. 5, we consider a scenario where 
all sparticles are too heavy to be produced at the LHC. A SU{5) test is possible through 
the SUSY-induced flavour-changing dipole operators. Beyond standard LHC searches, we 
propose a precision search for dipole operators in ultraperipheral collisions. As an aside it 
is pointed out that EW fluxes in p-p, p-Pb, Pb-Pb collisions are expected to be coherent, 
and thus drastically enhanced. 

In Sec. 6 we consider the natural SUSY case. Depending on the stop/wino ordering, a 
test on stop flavour-violation involving charm-tagging and a test involving top polarimetry 
are proposed. The latter is favoured in presence of a light mainly right-handed stop. In 
Sec. 7, both stops and scharms are assumed to be directly observable at the LHC. This top- 
charm SUSY scenario is found to be motivated in various 5d constructions. Tests relying 
on Higgs detection in cascade decays are presented. Summary and conclusions are given 
in Sec. 8. 

2 The sector of up-type squarks in SU{5) theories 

In this Section, we review the sector of up-type squarks in the light of the SU (5) hypothesis, 
and discuss a SU{5) relation which has, to our knowledge, only been reported in the earlier 
work of Ref. [33]. Starting from the superpotential given in Eq. (1.2), the lOjlOj term is 
symmetric, such that only the symmetric part of the Xu is selected. This enforces a 
symmetric top Yukawa coupling at the GUT scale, 

Uu = Uu- ( 2 - 1 ) 

Contrary to the aforementioned SU{5) relation of Eq. (1.1) which is potentially modified 
by model-dependent GUT threshold corrections or higher dimensional operators, Eq. (2.1) 
remains exact whatever extra corrections are included. More precisely, for any operator of 
the form OijlOiWj, the antisymmetric part of Oij always vanishes in the contraction with 
lOjlOj. Let us emphasize that yu is symmetric, but generally not hermitian. 

For a non-SUSY theory, Eq. (2.1) does not seem particularly exploitable. In this case 
the only physical parameters are mass eigenvalues and CKM angles, and this does not seem 
to be enough to find out whether or not yu is symmetric. This situation changes once we 
consider SUSY. More precisely, thanks to the required SUSY breaking, more degrees of 
freedom of the Yukawa matrices become physical. This is due to the presence of the scalar 
couplings, which in general are not flavour singlets. Just like the SU{5) relation Eq. (1.6) 
for the SUSY-breaking terms comes along with Eq. (1.1) in the sector of slepton and down 
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squarks, we have a relation between the trilinear SUSY-breaking terms in the up-squark 
sector associated to Eq. (2.1), 

an = On • (2.2) 

As we will see in the following, this property provides a potential way to test the SU{5) 
GUT hypothesis. We remind that we work in the assumption that the SUSY breaking 
source does not break SU{5), i.e. the SUSY breaking is parametrized by a spurion E-term 
which is singlet under SU{5). Going beyond this simplest possibility, one observes that 
a spurion E-term in the 24 and 75 of SU{5) would lead respectively to and 

ttu = a\j^. This comes from the tensor product 10 x 10 = Sgym -|- 45antisym + SOgym- We do 
not consider these possibilities here. 

A crucial point is that Eq. (2.2) survives when rotating the supermultiplets into the 
Yukawa matrices eigenvalue basis (generally referred to as the super-GKM basis). Gener¬ 
ally, the singular-value decomposition of the complex Yukawa matrices is yf = UjyjkEj, 
where Uf and Wf are unitary matrices and y/ is the diagonal matrix containing the (posi¬ 
tive) eigenvalues. The rotated trilinear terms are therefore aj = WjufUf. The fact that the 
Yukawa coupling is symmetric implies U* = Wu, such that we can write = W^yuWu- 
Further, we obtain = WuCLuW*, which is a symmetric matrix if Uu is. One can re¬ 
mark that both properties Eqs. (2.1) and (2.2) are necessary to keep Uu symmetric in the 
super-GKM basis 

Let us emphasize that the relations (2.1), (2.2) of our interest are confined to the 
flavour space of up-(s)quarks, and are therefore expected to be more stable against quantum 
corrections than the quark-lepton relations (1.6) and (1.7). Moreover, let us note that there 
is no connection between the SU{5) relation and potentially dangerous FCNCs 

in the down-sector. In our later analyses of SUSY scenarios, we assume some up squarks 
to be observed, which means that constraints from down-sector FCNCs are assumed to be 
fulfilled by construction. 

What we have found up to now is summarized by the implication 

|5C/(5)-type SUSY GUt} ^ |y„ = y(,, Uu = at the GUT scale} . (2.3) 

The next step is to see how this translates into physical properties, which can be tested 
at the TeV scale. Since the property given in Eq. (2.3) is confined to up-type squarks, we 
only need to scrutinize this sector. The object comprising all necessary information for our 
purpose is the up-squark mass matrix contained in the Lagrangian 

£ D -u^MIu, (2.4) 

where u = {ul, cl, ii, ur, cr, IrY contains the six up-type squarks. The 6x6 mass matrix 

is hermitian by construction. 

®More generally, Uu remains symmetric for any rotation of pu of the form j// = UjpfUf, i.e. preserving 
its symmetry. 
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In the broken electroweak phase and within the super-CKM basis, the up-squark mass 
matrix takes the form 


2_^ -^Q + 0(t>^)l3 ^ (1 + + ...) +0(n„MsusY)l3 

(1 + Cq,^ + ...) + 0(t’ti-MsusY)l3 + 0{v‘^)ls 

[ , ( 2 . 5 ) 

Here, the parameters Mq and are the scalar masses of the left- and right-handed up- 
squarks in the super-CKM basis, and Msusy is the typical SUSY scale. The Hu Higgs field 
of the MSSM is decomposed into its vacuum expectation value Vu and its light component 
h (the SM-like Higgs boson) with fraction Cq. The heavy neutral component H is not 
shown. More detailed expressions can be found in any SUSY review (see, e.g.. Ref. [37]) 
and are not needed for the rest of this work. 

We can see that the property clu = implies that the off-diagonal block of is 
symmetric. This property will be at the center of our attention in this work. Finding 
evidence for a symmetric off-diagonal block of would constitute a rather striking hint 
in favour of the SU{5) hypothesis. Finding evidence against a symmetric off-diagonal block 
in Yf| would prove the SU{5) hypothesis to be wrong. From this matrix we see that the 
physics needed to test this SU{5) relation will have potentially to do with up-squark flavour 
violation, chiral (LR) symmetry, and will possibly involve the Higgs boson. 

While the property is exactly valid at the GUT scale, it could be lost in the 

renormalization group (RG) evolution and not observable at C>(TeV) scales, just as it is 
the case for the SU{5) relations Eqs. (1.6) and (1.7). In the following, we are going to 
quantify how much information about the property au = disappears through the RG 
flow between the GUT and the TeV scales. 

Let us first consider the relevant RG equations [37]. We have to study the beta 
functions of the Yukawa and trilinear matrices, yu and a^, since these beta functions 
conserve the property au = any scale if they are symmetric. The one-loop MSSM 

beta function of the up-type Yukawa matrix is given by 



= Vu 


3 Tr{yi?/„} -h 3 ylvu + y\yd - y “ 3S'! - ^ 5? 


( 2 . 6 ) 


In this expression, the gauge and the trace terms are flavour singlets, and yuVuyu is sym¬ 
metric by hypothesis. The only non-symmetric contribution is the term yuUdVd- SUSY, 
the relative magnitude of the up- and down-type Yukawa couplings depends on tan P = 
such that the contribution of this term to the asymmetry grows with tan/3. Given the 
observed low-energy quark masses, one expects y^yd to compete with yhyu only at very 
high tan/3. Moreover, the yuV^Ud contribution is suppressed by the elements of the GKM 
matrix, as yd and yu would commute if they were simultaneously diagonalizable. One can 
thus expect that yu stays symmetric to a very good approximation at low energy. 
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The one-loop MSSM beta-function of reads 


167r^^a„ = ay 


3 Tvlylyu} -h 5 ylyu + y\yd - y 53 “ 352 “ ^ ffi 


+ yu 


32 26 

GTrjauy^} -|- 4,y\au + -|- — M| -|- 6M| -|- — Mf 


(2.7) 


Again, the gauge and trace terms are flavour singlets. But the other terms are generally 
not symmetric, because generally does not commute with ytyu, nor does ad with yuy\- 
Thus, au stays symmetric to a good precision at low energy only if the RG flow is dominated 
by gauge contributions. Therefore having a symmetric at low energy seems not to be 
such a generic feature. In practice, however, the beta-function is often dominated by the 
large gluino mass contribution. Moreover, as the M| contribution is positive, it decreases 
au with the energy, such that the non-symmetric terms become smaller and smaller in the 
beta function. Therefore, one can expect that although is generally non-symmetric, its 
asymmetry at the TeV scale remains rather small in many concrete cases (see Sec. 3). 

Overall, we expect and a„ to remain symmetric to a good approximation at scales 
that can be tested, e.g., at the Large Hadron Collider (LHC). Somehow, this S'?7(5)-type 
relation persists because it occurs inside the flavour space of the up-(s)quark sector, while 
the other 5C/(5)-type relations Eqs. (1.1) and (1.6) are relating different sectors, and are 
thus more sensitive to quantum corrections. 

The theoretical, model-dependent deviations from induced by the RG flow 

must be considered as an irreducible uncertainty. Quantifying this uncertainty by the 
implication Eq. (2.3) becomes therefore 

{517(5)-type SUSY GUt} ^ |d„ = a^ (l + 0(e”)) at the TeV scale} . (2.8) 


We numerically estimate in Sec. 3 for typical scenarios within the MSSM. 

Beyond the MSSM, one can also expect that yu and be symmetric at low-energy to 
a good precision. The condition is that the non-symmetric terms in Eqs. (2.6) and (2.7) do 
not dominate, and that the possible hidden sector contributions do not contribute either 
to For example, a hidden SUSY breaking sector would need to be flavour singlet to 
not spoil the symmetry. From now on, let us simply keep in mind that having <C 1 
is not absolutely generic, although it is present in many classic cases. In the rest of this 
paper we work in the super-CKM basis unless stated otherwise, and will not write the hat 
on au anymore. The low-energy SU{5) relation will be commonly denoted by ~ a^. 

Before coming to low-energy observables and to concrete numerical examples, let us 
set the notations of the up-squark parameters. The complete low-energy up-squark mass 
matrix in the super-GKM basis is denoted 


/m?i 


M- = 


^12 

m22 


•13 

mi4 


mi6 

2 

-23 

^24 

"^25 

^26 

2 

'33 

m|4 

"^35 



111.44 

mis 

ttiIq 




mle 






\ 


(2.9) 
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(860)2 

0 

0 

0 

(374)2 

0 

0 

0 

(412)2 



(1000)2 

0 

0 

0 

(3270)2 

(557)2 

0 

(557)2 

(1550)2 


{0‘u)ij 

0 

0 

0 

0 

0 

( 0^)23 

0 

( 0^)32 

-1500 


{o.d)ij 

0 

0 

0 

0 

0 

( 0^)23 

0 

{ad)z2 

400 


Mi^2 — 1000 

= (1000)2 

sign(^) = +1 


Table 1. Supersymmetry-breaking parameters at Q = Mgut of the MSSM reference scenario used 
in our numerical analysis. Masses and trilinear couplings are given in GeV. Two values for tan j3 
are considered, tan /3 = 10 and 40. 

We will denote the masses of the up squarks by rriui (* ~ with the convention 

niu^ < rriu^ < • • • < rriug. At the TeV scale, the squared mass matrix defined in Eq. 
(2.5) is diagonalized according to 

diag(m|j,...,m|J = R^MIrI, (2.10) 

where Ru is the 6x6 rotation matrix containing the flavour decomposition of each squark 
mass eigenstate Ui. 

3 Persistence of the SU{5) relation at low-energy 

In the previous Section, we saw that the SU (5) relation ~ remains true at low scale 
up to a relative discrepancy measured by the parameter This discrepancy is induced 
by the RG flow, and a qualitative study of the MSSM RG equations suggests that e*’’'’ 
is parametrically small in most of the parameter space. In this Section we numerically 
evaluate the magnitude of It is mandatory to know this magnitude as it constitutes an 
irreducible bound on the precision at which one can potentially test ~ A too large 
would imply that this relation is spoiled by large model-dependent quantum corrections. 
Here we present the results for a reference scenario of the MSSM. The scenario is 
defined at the GUT scale by the parameter values given in Table 1. The SUSY-breaking 
parameters of Eqs. (1.4) and (1.5) satisfy M^q = Mg = = Mj,, M| = = M|, and 

CLd = CLi Sit Q = Mgut- The only CP-violating phase is the one of the CKM matrix. 

We compute the SUSY spectrum in the DR scheme at the scale Q = 1 TeV according to 
the SPA convention [38] using two-loop RGEs. We use the spectrum calculator SPhenoMSSM 
obtained from SARAH [39] and SPheno [40]. To quantify the discrepancy between and 
a(,, it is convenient to use a quantity normalized with respect to the SUSY scale. Therefore 
instead of using we introduce the asymmetry normalized to tr(Al|), which is arguably 
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AMbs 

(17.76 ± 3.37) ps-i 


(1T0.25) 

BR(i30 ^ 

(3.12 ±2.08) X 10"^ 

BR(6 —)• 57) 

(355 ± 68) X 10"® 

BR(r —5- /i7) 

< 4.5 X 10"® 


Table 2. Low-energy flavour data taken into account in the numerical analysis, 
representative of the SUSY scale. More precisely, we define 

^ _ \{o,u)ij 

which is related to by Aij = |(a„)ij|/Tr{Yf?}^/^. Note that the asymmetry Aij 

being dimensionless, it is independent of the overall scale chosen for our reference scenario. 

There are in principle three asymmetries to be investigated; Au, A13, and M 23 . How¬ 
ever, the 1 — 2 and 1 — 3 mixing terms appear to be too constrained by the requirements 
that no tachyons should be present in the low-scale spectrum and that loop corrections to 
the SM Yukawas should not be too large. Moreover, the largest asymmetry is expected to 
come from the 2 — 3 sector, so that M 12 , -4.13 < - 423 . We focus therefore on the asymmetry 
-423- 

As discussed in Sec. 2, the main dependence of A23 is with respect to the parameters 
(«u)23,32 and (0^)23,32- We vary these parameters assuming that the SU { 5 ) hypothesis is 
true, i.e. ( 0^)23 = (a «)32 at Q = Mgut- We also take ( 0^)23 = ( 0^)32 for simplicity. 

For a complete analysis of the reference scenario, we take into account constraints from 
flavour physics at low-energy, among which the Bg mesons oscillation frequency AM^s 
and the ratio of the AT-meson mixing parameter normalized to the Standard Model value 
The list of constraints is given in Table 2. For each constraint, experimental and 
theoretical errors are combined in quadrature. The central values and experimental errors 
can be found in Refs. [41], while the theoretical errors are taken from Refs. [42-44]. In 
our reference scenario, the SM Higgs mass is found to be about nif^o = 124 GeV, which 
is consistent with the latest ATLAS result, nif^o = 125.36 ± 0.41 (sys-|-stat) GeV [45], 
combined in quadrature with a theoretical error estimate of ~ 3 GeV. 

We evaluate the asymmetry A 23 at the TeV scale for a scan in the parameters ( 0^)23 = 
( 0^)32 and (arf )23 = ( 0^)32 at the GUT scale. The numerical results are presented in Figs. 
1 and 2 for tan /I = 10 and tan /3 = 40, respectively. The thick black isolines indicate the 
asymmetry -423 at the scale Q = 1 TeV in percent. Note that the bounds | ( 0 ^) 23 1 < 600 
GeV and [( 0 ^) 23 ! < 400 GeV are taken into account in order to avoid tachyons in the 
resulting mass spectrum. The leading constraints are BR(i3° —)■ /r^), AMbs (see Fig. 1), 
and BR (6 —> 57 ) (at high tan/3, see Fig. 2). We find e;y/e|^ ~ 1 due to the fact that 
we have chosen real SUSY-breaking parameters. We have numerically verified that this 
observable constitutes a stringent constraint on CP-violating phases. Finally, let us note 
that for both low and high tan /3, BR(r —)• //y) never exceeds the upper experimental value. 


i A j , 


(3.1) 


0=1 TfiV 
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(a'u)23,32(0 = Mqut) 


Figure 1. The asymmetry ^ 2,3 (black solid line) together with the 2cr exclusion bands from AMbs 
( blue dashed lines) and BR(_B° —)• fifi) (red dashed-dotted lines) evaluated in the reference scenario 
of Tab. 2 for the case of low tan/3 = 10. The grey area represents the allowed zone once all 
constraints are taken into account. 

Moreover, the value of the Higgs mass also does not strongly depend on the variation of 
(a«)23,32 or (0^)23,32- 

Coming back to our main center of interest, we can see in Figs. 1 and 2 that the 
asymmetry ^23 does not exceed a few percent. As expected from RGE considerations, we 
can conclude that none of the asymmetries Aij exceeds this value on a large part of the 
MSSM parameter space. In the following, we find that, during the LHC era, such a level 
of precision will be most probably difficult to reach. We will therefore not mention or 
Aij for the rest of this Paper. 

4 Strategies and tools for testing the SU{5) hypothesis at the TeV scale 

Any strategy that can be set up to test the SU{5) relation necessarily relies on a 

comparison involving at least two up-squarks. Apart from this relation, the squark matrix 
is in general arbitrary, so that each of the six up-type squarks can take any mass. Some of 
the squarks can be light enough to be produced on-shell at the LHC, while others may be 
too heavy such that they appear only virtually in intermediate processes. 

As a result, a panel of possibilities for SU{5) tests will appear, depending on the 
exact features of the up-squark spectrum. It can be convenient to split the possibilities 
of SU{5) tests into three categories, depending on whether the tests involve only virtual, 
both virtual and real, or only real up-type squarks. In the following, we outline the tools 
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H 





M23m{Q = Mgut) 


Figure 2. Same as Fig. 1 for the case of high tan/3 = 40. The green dotted line represents the 2cr 
lower bound on BR (6 —>■ 57 ). 


that are relevant for each of these cases. Note that the SU{5) tests on virtual up-type 
squarks will necessarily involve one-loop processes, because the interaction terms probed 
in the tests are i?-parity conserving. In contrast, the SU{5) tests on real squarks can in 
principle rely on tree-level processes only. 

Dealing with the eigenvalues and rotation matrices of the full 6 x 6 mass matrix is 
in general rather technical, and may constitute an obstacle on the quest for simple SU{5) 
relations among observables. However, depending on the pattern of Ad|, two expansions 
can be used in order to simplify the problem: the mass insertion approximation (MIA) and 
the effective field theory (EFT) expansion. 

The MIA is an expansion with respect to small eigenvalues-splittings, and is in par¬ 
ticular valid in presence of small off-diagonal elements. In contrast, the EFT expansion 
can be applied if a large hierarchy is present in Ad?. The heavy up-type squarks with 
mass matrix are integrated out, the expansion parameter is and can have 

arbitrary off-diagonal entries. The two expansions, which we detail in the following, are 
therefore complementary. 

It is clear that the feasibility of the SU (5) tests we will setup depends crucially on the 
amount of data available - whatever they involve real or virtual up-squarks. This feasibility 
needs to be quantified using appropriate statistical tools. Whenever a 517(5) test can be 
obtained through a definite relation among observables, we will use a frequentist p-value 
in order to evaluate to which precision this relation can be tested, for a given significance 
and amount of data. We will systematically report the expected precision associated with 
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this test. 


4.1 Tree-level up-squark effective theory 

The pattern of the squark masses can in principle exhibit some hierarchy, and this actually 
occurs in many classes of new physics models. In such a situation, the physics of the light 
squarks can be conveniently captured into a low-energy effective Lagrangian, where heavy 
squarks are integrated out. This has, e.g., be discussed in Ref. [33] for the situation where 
only two squarks are light and the remaining four can be integrated out. The higher¬ 
dimensional operators present in this tree-level Lagrangian will serve as a basis for the 
SU (5) tests based on both virtual and real squarks. 

Following Ref. [33], the up-squark mass term of Eq. (2.4) can be reorganized such that 

£ 3 n'Mlu ^ , (4.1) 


where cj) contains the heavy states and cj) the light ones. The relevant piece of the corre¬ 
sponding Lagrangian has the general form 

C D + + , (4.2) 


where the operators O and O represent the interactions with other fields, that are poten¬ 
tially exploited to probe the up-squark sector. 

Assuming that the matrix norm of is large with respect to the energy at which 
the theory is probed, i.e. ||M^|| = (X^j^ ^ the heavy squarks (j) can be 

integrated out. We are thus left with the low-energy effective Lagrangian of light squarks. 


I -k h.c. . (4.3) 


-k 


O 


O - 0(M"2 _ 


We keep only the leading and the subleading terms of the expansion relevant 

to build SU{5) tests. This effective Lagrangian contains in principle higher dimensional 
couplings and derivative terms, which are either subleading or irrelevant for the observables 
we are going to consider, and are thus neglected. Eq. (4.3) has been obtained using the 
field redefinition —)■ (1 — in order to canonically normalize the light 

squarks kinetic terms. The { ,} denotes the anti-commutator. 

Erom Eq. (4.3) we see that flavour-violating couplings of the light squarks enter at 
first order and are controled by The flavour-conserving couplings will instead 

be modified at the second order. The mass matrix M^, associated to the light squarks, 
receives a correction independent of at first order, and corrections proportional to 
at second order. The imprint of the heavy up-squarks in the Lagrangian of Eq. (4.3) 
appears as corrections to the masses and couplings of the light up-squark states. Physically, 
these corrections have to be understood both as tree-level exchange of heavy up-squarks. 


- 12 - 






and as the first terms of the expansion with respect to the small parameters that describe 
mixing of heavy and light squarks. While in Ref. [33] it has been applied to the special 
case of two light squarks, this framework is more general and can be used for any number 
of states that are accessible at the LHC, while the others are integrated out. Finally, we 
remind that, although this effective theory approach might bear some resemblance with 
the mass insertion approximation, the two approaches are fundamentally different. 


4.2 One-loop effective operators 

When all the up-squarks are heavier than the typical LHC scale, they can only appear 
off-shell. Potential tests of the SU{5) hypothesis are thus based on virtual squarks only. 
All squarks can be integrated out, and the resulting low-energy Lagrangian then contains 
the SM plus possible other light SUSY particles. Using other light SUSY particles to 
test Qu ~ does not seem to provide attractive possibilities. In particular gauginos are 
flavor-blind, and bottom squarks depend on their own flavour structure. 

Here, we focus instead on the case where the whole SUSY spectrum is heavy. The 
resulting low-energy Lagrangian then has the form 

CeS = CsM + + 0(1/M|usy) • (4.4) 


The effective operators with dimension larger than 6 are not relevant for our present study. 
The dimension-six Lagrangian takes the form 


£( 6 ) 


E 


a,: 


^'^SUSY 


Oi, 


(4.5) 


where the are dimensionless constants accompanying the operators Oi. This Lagrangian 
contains 59 operators, which have been first fully classified in Ref. [47]. 

In SUSY models with conserved i?-parity, all CP-even operators are generated at loop 
level. In general, they receive contributions from a large number of loop diagrams, and we 
shall identify which among these operators can be potentially relevant for a SU{5) test. 
We are thus looking for one-loop diagrams involving 


• at least one {au)ijhuqiUj vertex, as it encloses the SU{5) relation Eq. (2.2) that we 
want to test. 


• fermions on the outgoing legs in order to access information on the flavour structure 
of the {au)ij coupling. 


These requirements naturally lead to the sector of up-type flavour-changing dipole opera¬ 
tors of the form 


O. 


V 

u,ij 


= qia^^u 


'j Hu 




(4.6) 


where V = {G,W,B) and = |[7^,7^]. The key point to obtain a test of the SU{5) 
hypothesis, is to be able to distinguish between and The generation of these 

operators and the possibilities they offer to build a SU{5) test will be discussed in Sec. 5. 
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4.3 Mass-insertion approximation 


If the relative difference between the eigenvalues of the mass matrix Al| is small, an 
expansion can be made by taking this difference as the expansion parameter (see, e.g.. 
Refs. [48-50]). This is the so-called mass-insertion approximation (MIA) In the super- 
CKM basis, the mass matrix is written in terms of the mass insertions du^ {ct, = L,R) 
as follows. 




m? 


1 + 




(4.7) 


where is an average squark mass, that can, e.g., conveniently be chosen to be = 
Tr{A4|}. Note that, here, the <5“^ are 3x3 matrices in flavour space, and we have 
since the mass matrix Al| is hermitian. The SU{5) relation ^ 
translates on the chirality-flipping, flavour-violating mass insertion matrices as 


^LR 

'-'u 




(4.8) 


or equivalently, in terms of the generation-mixing entries, 

= (<5«'^)i3, {^u^) 32 = {Su^)23- (4.9) 

The MIA also applies when only a subset of the Al| eigenstates is nearly-degenerate. 

The MIA is commonly used in the broken electroweak phase with (H) = vj^/2 ^ 0. 
However, the effect of electroweak breaking in the mass matrix is small by assumption, as it 
participates in splitting the eigenvalues. This is in particular verified whenever MgusY 
u/\/2. Therefore electroweak breaking does not have influence on the mean squark mass 
rriq by assumption. As a result, the MIA can also be used in the unbroken electroweak 
phase. In this case, the VEVs have to be replaced by the original Higgs fields. Electroweak 
symmetry breaking appears in the squark mass matrices exclusively through the chirality- 
flippling mass insertions 6 ^^, 6^^. More precisely, the latter are proportional to the VEV 
according to y/Msusy, while the chirality-conserving ones are not. Thus, the use 

of the parameters 5^^ and 6 ^^ remains unchanged. In contrast, for the chirality-flipping, 
flavour-violating mass insertions, which we are interested in, we have to replace 

^ ( 4 . 10 ) 

Vu 

for i ^ j. This is also a way to recover the physical Higgs boson in the complete mass 
matrix Eq. (2.5). 

4.4 Expected precision 

When the SU{5) test consists of a simple relation, for example among event rates at the 
LHC, one can use a frequentist p-value to evaluate the potential of the test. Assume that 

^Note that the MIA is sometimes considered as an expansion with respect to small off-diagonal elements 
of At|. This is a more restrictive expansion, which is in general not equivalent to the MIA as we define it 
here [49, 50]. 
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there is a relation R{Oi) among observables quantities Oi, that satisfies 

R[R{Oi)\ = 0 (4.11) 

whenever the SU{5) hypothesis is verified. Here, E denotes the expectation operator. 
We use a p-value test to quantify whether or not the observed value of R, denoted f, is 
compatible with zero. This test is not trivial to perform in general, but simplifies whenever 
the probability density function (PDF) of f can be approximated by a half-normal law 
f ~ 2M{r,a)6{r) (see e.g. Ref. [51]), where 9 is the unit step function. 

In this case, the compatibility of r with the r = 0 hypothesis can be expressed in terms 
of a statistical significance as 

Z = ■^. (4.12) 

a 

Here, Z = 4>“^(1 — p) translates the p-value into a significance given in terms of standard 
Gaussian deviations. is the standard Gaussian repartition function. 

For a given relation satisfying Eq. (4.11) and an hypothesized amount of data, one can 
evaluate the expected value of a. For a given significance Z, the quantity 

Pz = Za (4.13) 

represents the expected precision to which the relation Eq. (4.11) can be tested. Throughout 
this paper we will systematically report the value of Pz when a relation like Eq. (4.11) is 
available. 

As an example, consider the relation aXi = 6 A 2 . The Xi are assumed to be normal 
variables, centered on pn and with variance af. Defining R = \aXi — 6 A 2 I, the PDF of 
R takes the form //j ~ 2J\f{r,a‘^)6{r) where r = api — bp 2 , and the variance is given by 
(T^ = a'^af + b‘^o' 2 - In order to compute the expected significance defined in Eq. (4.13), 
we use this value of a assuming that the SU (5) hypothesis is true. Note that this implies 
r = 0 by construction. 

In practice it is useful to normalize R such that E[i2] G [0,1]. The expected precision 
can then be expressed in percent. Concretely, having Pz = 20% means that a violation of 
the relation E[i?] = 0 by 20% or more (i.e. E[R] > 20%) can be assessed with a significance 
of Za. The use of a normalized quantity is more convenient to interface the SU{5) test 
with other studies - either experimental or theoretical. Normalizing R implies that one 
frequently encounters the random variable \A — B\/{A + B), where A and B are potentially 
correlated. We will use the notation E[A] = pA, E[-B] = pb, V[A] = a\, V[H] = cr^, 
Cov[A, B] = pa AO’s- The following relations are useful in order to derive the expected 
precision, 

P[A/B] 

X[A/B] 

X[{A-C)/{A + C)\ 


E[A]/E[B ], 


Pa<7b + Pb^^a - ^paACTBtJ-AtJ-B 


^ p\<tI + pI<t\ 

{pA + Pc)* 


(4.14) 

(4.15) 


4V[(A)/(A + C)] with Cov[A, C] = 0 . 
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These approximations are in particular valid in the Gaussian limit, where a ratio of random 
variables is described by the Fieller-Hinkley distribution [52], For a large number of events 
Na = HA ~ cta similarly for C, one has for example 


Y[\Na-Nc\/{Na + Nc)] 


4NaNc 
{Na + Nc)^ • 


(4.17) 


5 Heavy supersymmetry 


We can now apply the strategy and tools presented above in various SUSY scenarios. In 
this Section we assume that the SUSY scale is large with respect to the scales probed by the 
LHC. This of course increases the fine-tuning of the electroweak scale which is a primary 
motivation for low-energy SUSY In such a situation it is appropriate to integrate the 
whole SUSY sector including the charged Higgses. As discussed in Sec. 4, the dimension-six 
operators that one should scrutinize to test the SU{5) hypothesis are the up-sector dipoles 


/’( 6 ) 

■'^l-loop 


a. 


D 




^'^SUSY 




a, 


w 

V 


/\/f2 

^'^SUSY 




a. 


/Vf 2 

^“^SUSY 


qiC^ UjHuB^i/ , 


(5.1) 

where i ^ j. Our ability to build a S'[/(5)-test relies crucially on distinguishing the chiral¬ 
ities (i.e. between L and R) at some level. The only quark for which this can be done is 
the top quark, as the lighter quarks hadronize too fast. Our test has therefore to rely on 
top polarimetry, which itself relies on the V — A structure of the leading top decay t —)• bW 
(see e.g. Ref. [54]). 

We therefore focus only on the dipole operators involving a top quark, 


AG) 

-^l-loop 


D 


a 


^'^SUSY 




a. 


fJ.i' ' 

^'^SUSY 


ULia^^''THRH^Gl 


{lU 






a 


w 

3i 


^ ^SUSY 


tLCT^’^URiH^Wj^, + 


Oi 


w 

i3 


IJ-i' ' 2 /Vf2 

^ ^'^SUSY 


UL^a^•^ThRHuW, 


3 


a 


^'^SUSY 


tia^'^URiHuBau + 


a. 


' m 2 

^'^SUSY 


ULi(T^''tRHuB^ 


flly ? 


(5.2) 


where i = 1,2. 

The operators given in Eq. (5.1) are generated at one-loop by penguin diagrams. They 
receive their main contributions from the chargino-down-squark, charged Higgs-down- 
quark and gluino-up-squark amplitudes. The latter contribution is expected to dominate 
because it is enhanced by the QCD coupling. Moreover, contrary to the well-known case of 
down quarks, higher-loop corrections in the up-sector from flavour-changing self-energies 
[55, 56] are tan/3-suppressed, so that they cannot become large. We therefore end up with 
dipole operators mainly generated by the gluino loop. 

®On the other hand the SUSY scale can still be far from the Planck scale, so that it still improves a lot 
the gauge-hierarchy problem. To go further into naturalness and hne-tuning arguments, these notions need 
to be properly defined and quantified, see Ref. [53]. 
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From now on we assume that the up-squark mass matrix is degenerate enough so that 
the MIA applies. The expression of the electroweak and chromo-dipoles operators are then 
given by 



9s 

(^V2 

mg6i^FP{xg] 

^'^SUSY 

167r2 m? 

Vu 240 


55s / 

^ V2 - 


A//'2 

^‘^SUSY 

IbTT^m? \ 

Vu 30 


g'9 s ( 

V2'- 

^g^ki^^EW ) 

^‘^SUSY 

IbTT^ m? ( 

Vu 180 


- 5 


+ 


+ 


yt6g^FP{x,) 
36 

9 ) 

54 i 


(5.3) 

(5.4) 

(5.5) 


and the a^s coefficients are obtained by replacing 5^1^ —)> 6 ^/" —>■ 6 ^^. Here, we have 

defined Xg = rrig/mg, the form factors satisfy F{1) = 1, F{x) —)• 0 for x —)> oo, and are 
given in App. A. Moreover, niq denotes the average squark mass introduced in Sec. 4.3. 
The above expressions have been deduced using loop functions from Refs. [57, 58] and by 
making use of Sec. 4. The coefficients of the EW operators in the broken phase. 


'y Z 

4®^ioop ^ -vS^iL(r>^"^RiHuF^, + -^^iLa>^''uRiHuZ^,, (5.6) 

^'^SUSY ^^SUSY 

are obtained by replacing aj^ —)• /2 -|- Cwafj and —)■ /2 — SwO-fj- Here Su, and 

Cyj are the sine and cosine of the weak angle. 

As already pointed out above, for the purpose of obtaining a SU{F) test, the contribu¬ 
tions that we are interested in are the chirality-flipping ones, 5^^. The chirality-conserving 
contributions from 5^^ and 5^^ are suppressed by a factor ytVu/mg = mtsp/rrig) with 

respect to the contributions from 5^^. This factor is compensated if one considers that 
^LR 

oc vau/rn^, as it is usually the case in SUSY models with soft breaking. However, 
the magnitude of the a„-term can be larger than mq (for example from naive dimensional 
analysis [59]). One should therefore let 5^^ rsj 0(1). 

Considering usual models, one should also notice that 



(5.7) 


at low-energy to a very good aproximation. We can now identify the consequences of our 
517(5) relation ps at the LHC. This relation implies the following equalities between 
the coefficients of the dipoles in Eq. (5.2), 


aS 




w 

“3* 


W 

OiiS : 




ai- 


(5.8) 


Concretely, the dipole operators induce both new top decays and top production modes 
at the LHC. These processes will be discussed in the Subsecs. 5.2 and 5.3 of the present 
Paper. 
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5.1 SU{5) test through single top polarimetry 

We have seen above that measuring the top spin is a necessary ingredient to build a SU{5) 
test that distinguishes between the and dipole operators. The expected precision 
associated with such polarization-based tests can be evaluated in a generic way. The top 
spin has to be measured with some polarization-sensitive observable z with distribution fz- 
For the t —)• bW decay, z can be for example the angle between the top and the lepton from 
the IT-decay, or the 5-quark energy in case of a boosted top [54]. The usual way to proceed 
to get information on the top spin is to split the phase space into two subsets V = +V^. 

The total sample of N events is then split into two subsamples N± = N\z£t>± satisfying 


N = N+ + N_, E[Ar±] = E[iV] [ fziz)dz. 

Jv± 


(5.9) 


The choice of 'D± is in general a freedom of the analysis, although it is preferable to have 
N+ rsj N- to minimize the statistical error. In what follows, we will systematically choose 
the subsets of phase space V± such that E[A^+] = E[A^_] if the SU (5) hypothesis is satished. 

The power of the SU{5) test depends on the power of the spin-analyzing variable, 
which is parametrized by k G [0,1]. Here we assume that the z-dependent distribution 
used to analyze the spin has the form fzoc{l + KPtz), where z G [—1,1] and Pt = ±1 
is the top polarization [54]. A similar approach can be done for any other distributions. 
The lepton angle in t —)• bW has a maximal spin-analyzing power n = 1 because of the 
V — A structure of electroweak interactions. In that case 2 ; is the top-lepton angle cosine 
in the top rest frame. Note that n < 1 from unitarity. Eor the subdomains of z we choose 
P_i_ = [0,1], P- = [—1,0[. Note that this matches the usual definition of a forward- 
backward asymmetry. 

In order to quantify the precision at which the SU (5) relation can be tested, we intro¬ 
duce the normalized asymmetry 


which satisfies 


R = 

2 \N+ - N-\ 

K N+ + N- ’ 

(5.10) 

E[i2] = 

llasiP - |ai3p| 

(5.11) 

asiP + ajsP 


Assuming B background events for the total signal, the expected precision defined in Sec. 
4 reads 


K JM 


(5.12) 


Note that Pz 00 for k —)• 0, as no spin information is available in this limit. The 
expected precision is shown in Eig. 3 as a function of the background and total event 
numbers B and N. The SU{5) hypothesis starts to be testable when Pz is lower than 


100 %. 

Eor K = 1 and B = 0, we see that the hypothesis starts to be testable for > 16 at 
2(7 and for A" > 36 at 3a significance. About 144, 900, and 3600 events are needed to test 


the relation at respectively 50%, 20%, 10% precision with 3a significance. Clearly, and 
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Figure 3. Expected precision Pz for top-polarization based SU{5)-tests as a function of the total 
number of signal and background events N and B. The spin-analyzing power is fixed to k = 1. 
Plain (dotted) lines denote 2cr (3cr) significance, respectively. Blue, purple, orange, and red lines 
respectively show ^ 2,3 = 100%, 50%, 20%, 10% isolines of expected precision. 


contrary to a mere signal discovery, a substantial amount of signal events is necessary in 
order to test the SU (5) hypothesis. 

5.2 Existing LHC searches 

Let us now discuss the various LHC processes that we expect to be relevant to perform the 
SU{5) test described previously. The dipole-induced anomalous top couplings induce the 
two-body decays 

t^qZ, t^qg. (5.13) 

Three and four-body decays also exist, they are obtained from the ones above by adding an 
extra Higgs in the final state and replacing Z by W~^W~ Among the two-body decays, 
t —)• qZ and t ^ qj carry information about the dipoles, that can be accessed through 
polarization of the outgoing gauge bosons. The t —>• (77 process has been searched for at 
the Tevatron and LHC [60, 61], the leading CMS bound is provided in Table 3. 

We now turn to dipole-induced top production. All the LHC processes we can consider 
have the particularity of featuring a single top. For a proton-proton collider, the main 
partonic processes one can think about are g q ^ t, g q ^ tZ, gq ^ t'y, gg ^ tq. Apart 
from the latter, various LHC searches have already been performed in these channels, see 
[62-66]. The sensitivities translated on dipole coefficients are given in Tab. 3. The leading 

®As the dipole operators of Eq. (5.2) involve a Higgs field, the effective vertices that they produce can be 
either taken with (H) = vj'/i, i.e. the tuV vertices, or with an extra Higgs boson, i.e. the tuVh vertices. 
The Higgs vertices can be overlooked if one directly introduces the effective operators in the broken EW 
phase, i.e. with {H) = vly/2, as customarily done in the anomalous-top coupling literature. 
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Dipole coefficients combinations 

95% CL limits 

SUSY 

( “31 P ^ -^SUSY 

< 0.19 TeV-2 [62] (CMS) 

7 • 10"® TeV-2 

(I«32P + I« 23 P)^^^^SUSY 

< 0.65 TeV-2 [62] (CMS) 

7 • 10"® TeV-2 

(l« 3 lP + l«f 3 p)^''^^SUSY 

< 0.68 TeV-2 [63] (CMS) 

1 • lO"'^ TeV-2 

( “32P + “23P) ^ ^SUSY 

< 3.44 TeV-2 [63] (CMS) 

1 • 10-“^ TeV-2 


< 0.029 TeV-2 [64] (ATLAS) 

3 • 10"^ TeV-2 


< 0.063 TeV-2 [64] (ATLAS) 

3 • 10"“^ TeV-2 


Table 3. Leading experimental limits on dipole operators. All the limits are given at 95% confidence 
level. The SUSY values are given for Msusy = 1 TeV, d ~ 1, and tan^ = 5. 

bounds on on^as, come from CMS and ATLAS searches for rare top decays, 

top -7 production and single top production. 

We observe that all these bounds seem to be far from the sensitivity required to observe 
the dipoles from SUSY, which are loop-generated. The typical order of magnitude for the 
SUSY dipole operators (Eq. (5.5)) indicated in the last 

column of Table 3. It might be in fact possible that none of the searches above will ever be 
able to reach the SUSY dipoles at the LHC, even with 3000 fb“^ of integrated luminosity. 

For illustration, let us consider the ATLAS gq ^ t search discussed in Ref. [64]. Know¬ 
ing the expected background B ss 10® and the 95% confidence level bound on |a:^p-|- 
for L = 14.2 fb“^ one can readily estimate the event rate by interpreting the significance 
as Z = S/\/B + S We obtain a « 52(|a^p -|- /M^^^y pb“^. Using the typical 

SUSY-induced values of shown above and extrapolating to L = 3000 fb“^, we find that 
only 28 signal events would be collected in this analysis. The background would need 
to be drastically reduced in order to obtain a 0(1) significance. A mere discovery of the 
SUSY dipoles using these analyses being apparently difficult, applying top polarimetry is 
even more compromised, as a substantial amount of events is necessary, as shown in Fig. 
3. 


5.3 Proposal for ultrapheripheral searches 

Given that the backgrounds are rather large in the LHC analysis described above, we now 
propose an alternative type of measurement relying on ultraperipheral collisions (UPCs). 
These are the quasi-diffractive processes that occur when the incoming protons or heavy- 
ions interact with an impact parameter larger than the sum of their radius (see Ref. [67] 
for a review). In the ultraperipheral regime the remnant of these nuclei are only slightly 
deflected and proceed in forward direction. The protons have a good chance to remain 
intact and can be tagged in the forward detectors planned for both ATLAS [ 68 ] and CMS 
[69]. Tagging both protons in fully exclusive processes provides a powerful background 
rejection, because the complete kinematics of the events are known (see e.g. Refs. [70-73]). 

^The ATLAS statistical analysis is much more evolved, but we expect this estimation to be enough for 
qualitative considerations. 
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Heavy-ion collisions have no pile-up, and do not need forward tagging. For both proton 
and heavy ions, UPCs are therefore an attractive tool for precision physics. 

As an accelerated charged object, the proton radiates short pulses of electromagnetic 
field. The associated photon flux emitted by the proton is proportional to the squared 
electric charge, oc ~ 0.1. Recent prospects for photon fusion in UPCs can be found in 
Ref. [74]. In the case of heavy ions, a striking feature is that all charges inside the nucleus 
radiate coherently, such that the photon flux scales as oc This phenomenon is 

known and exploited since a long time (see, e.g.. Ref. [75] and references therein). The 
enhancement factor is, e.g., 82^ ~ 4.5-10^ for Pb-Pb collisions. This enhancement has to be 
contrasted with the low expected heavy-ion run luminosity, typically 1 nb“^ per year. This 
type of experiment is useful if the large photon flux can compensate the small luminosity. 

We now point out that the same reasoning can qualitatively apply to electroweak 
charges. The W- and Z-boson fluxes can be treated similarly to the photon [76]. In 
the UPC regime, contrary to a seemingly widespread belief, electroweak charges radiate 
coherently just like the electric ones. Although the W and Z boson masses are responsible 
for the “short range” of the EW interaction in a non-relativistic view, in relativistic QFT 
these masses have to be put in balance with the centre-of-momentum energy of the di¬ 
boson collision. The maximum energy from a flux is estimated to be roughly j/R. For 
pp collisions, the maximum energy of the vector collision system is \/Wmax ~ TeV 
[74]. This is large with respect to mw or mz, so that the flux of transverse lU or Z is 
essentially identical to the one of the photon. Also, by the equivalence theorem, the flux 
of longitudinal W- and Z-bosons can be replaced by the flux of corresponding Goldstone 
bosons. For heavy-ion collisions, we have \/Wmax ~ GeV and ~ 160 GeV respectively 
for p-Pb and Pb-Pb, such that -v/Wmax ~ 0{mw,z)- In this case the EW fluxes are thus 
expected to be suppressed in comparison with the photon flux case. Understanding more 
precisely the features of these EW fluxes deserves a dedicated analysis. Note that with 
lead ions in initial states, the emission of W~^, W~ transforms Pb into Bi and Tl ions 
respectively. Such isotopes are unstable, but are long-lived with respect to the LHC scale. 

Here we simply focus on the effect of the total EW charge, i.e. the overall prefactor 
of the W- and Z-fluxes. For coherent fluxes f^r± and fz, we have /iy+ oc 2g^ ~ 0.9, 
flY- cc <7^/2 ~ 0.2, and fz oc /c^{l/4 — s‘^ + 2sf^) ~ 0.07. For the proton, the EW fluxes 
are therefore slightly enhanced for lU-bosons because of the magnitudes of the involved 
couplings. For the heavy-ion EW fluxes, in contrast to the electromagnetic case, the 
neutrons inside the nucleus are charged under the weak interaction. We obtain 


fw+ oc 
fw- oc 
fz oc 


g^z + Af 



2 

g^{2A-Z)^ 


2 


'{2Z -Af 


Z{2Z-A)sl + 2Zhi 


(5.14) 

(5.15) 

(5.16) 


For Pb-Pb collisions, we have therefore an enhancement of 4.1 • 10® for and 9.4 • 10® 
for with respect to the p-p case. For WW the enhancement factor is 200 times 
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Figure 4. Ultraperipheral processes potentially relevant to search for the dipole operators 
at high precision. The processes shown are the ones relevant for SU{5) tests. The vertex marked 
in red denotes the SUSY dipole. 

larger than the enhancement for 77 . Taking into account the gauge couplings, the factor 
from /^+/^y- is 10^ times larger than the factor from f^. Provided that the W fluxes are 
not strongly suppressed, processes involving intermediate lU-bosons in heavy-ion collisions 
at the LHC are thus rather attractive. 

The proton cannot emit a single gluon and stay intact because of color conservation. 
For UPCs however, an intriguing possibility is that the gluon comes from a Pomeron emit¬ 
ted by the proton [77] . In that case the proton remains intact and can be tagged. However, 
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the Pomeron itself cannot remain intact, such that the process is not fully exclusive. Taking 
all features into account, we now propose a selection of processes that deserves particular 
attention in order to probe the dipole operators with high sensitivity. The processes 

are pictured in Fig. 4: 

• "yj^tqin p-p collisions (Central exclusive production) 

The two protons can be tagged and the reaction is exclusive. However, because of 
the suppression for photon fluxes, the cross-section might be too small to detect 
the dipoles. 

• gg ^ tq through Pomeron fusion in p-p collisions (Double Pomeron exchange) 

The two intact protons can be tagged. However, this process is not fully exclusive as 
the Pomerons break up. A detailed study using diffractive PDFs would be necessary 
in order to evaluate the potential of this channel. This process is sensitive to the 
0^'"< dipoles. 

• Jg ^ tq in p-p collisions (Diffractive photoproduction) 

The proton radiating the photon can be tagged. One may also require a second 
intact proton. The remnant of the subsequent Pomeron would then be observed in 
one hemisphere only. This process is sensitive to , and may be an interesting 
trade-off between the two above processes. 

• WW tq in Pb-Pb collisions (Central exclusive production) 

This configuration fully relies on the huge enhancement factor from the coherent 
action of the EW charges of the nucleons. A computation of the IT-boson fluxes is 
required. This process is sensitive to the dipoles. 

• gi ^ g^ tq' in p-Pb collisions (Single diffraction dissociation) 

No proton tagging is necessary in this case, so that ones lets the gluon be emitted by 
the proton in order to take profit of the large gluon PDF. The luminosity is higher 
than in the Pb-Pb case, about 200 nb“^ per year. One takes profit of the 7 - or 
PF-flux enhancement from the heavy ion. This process is sensitive to the dipoles. 
Therefore this configuration might be a good compromise between p-p and Pb-Pb 
collisions. 

All these processes feature a top-quark plus a jet in the final state. Would the expert 
reader obtain an estimate for one of these processes, the potential to test the SU{5) hy¬ 
pothesis through top polarimetry could be directly drawn from Fig. 3. Identifying the jet 
charge - if possible - would be useful to reject fake QCD events [78]. Finally, we remark 
that all these processes are relevant for generic dipole operators searches beyond SUSY. 
Also, the enhancement of kF-boson fluxes from the total EW charge has implications for 
the searches of other anomalous couplings like triple and quartic gauge interactions. 

6 Natural supersymmetry 

We now investigate SU (5) tests in the framework of “natural” (or effective) supersymmetry. 
The mass spectrum of this scenario features a first and second generation of up-squarks that 
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are considerably heavier than the squarks of the third generation. This pattern is motivated 
by both naturalness considerations [79] and by ultraviolet (UV) constructions like partially- 
supersymmetric composite models [80] or supergravity contributions generically present in 
five-dimensional models [81]. 

Clearly, in natural SUSY, the up-squarks effective theory (see Sec. 4) will consist of 
two mostly stop-like squarks. Their mixing is not constrained and can potentially be large, 
which is often needed to satisfy the constraint coming from the observed Higgs-boson 
mass. The effective operators that appear when integrating out the heavy up-squarks can 
potentially induce flavour-changing stop decays, and violation of unitarity relations in the 
stop chiral decays [33]. 

All the information about the SU{5) relation ~ is enclosed in the higher¬ 
dimensional operators of the effective Lagrangian for the stops. Potential 517(5) tests 
have to rely on the effect of these effective operators in the stop sector. Therefore these 
tests involve both real squarks (the stops) and virtual squarks. We assume that both stops 
(U and ^ 2 ) are produced at the Large Hadron Collider (LHC). The total production cross- 
sections of stop pairs at next-to-next-to-leading order (NNLO) are = 89 fb“^ ± 10%, 
8.6 fb“^ ± 15%, 0.71 fb~^ ± 20%, and 0.08 ± 25% fb”^ for stop masses of = 700, 1000, 
1400, and 1800 GeV, respectively [82]. 

In the following, we develop 517(5) tests that already have been sketched in Ref. [33]. 
We refer to this paper for certain details of the derivation. We evaluate the potential 
of these tests using the frequentist treatment described in Sec. 4. We consider two cases 
concerning the mass ordering between the stops (m^^ ^), the bino (m^), and the wino (m^). 
We assume that the lightest neutralino is mostly bino-like, ~ H, and the second-lightest 
neutralino is mostly wino-like, X 2 ~ it is typically the case in scenarios with gaugino 

mass unification at the GUT scale. 

6.1 The ^ case 

This mass ordering allows the stops to decay either into the lightest neutralino ^ ^ B 
or into the second-lightest neutralino X 2 ~ order to build a 517(5) test, we are 

interested in the flavour-changing decays 

i —>■ W ufc —>■ BZjhulc and t —>■ Bufc. (6T) 

We assume that the stop masses are unknown, and that only the event rates of the stop 
decays into binos and winos, respectively denoted by Ny and Nl, are experimentally 
accessible. Moreover, we assume that a certain fraction, denoted by Ny and iV£, of these 
events can be charm-tagged. The remaining events Ny = Ny — Ny and = Nl — iV£ 
then contain both up-quark events and miss-tagged charm jets. 

Assuming the same charm-tagging efficiency Cc for Ny and N^ we have 

E[N^] = ecjyE[Ny], E[iV£] = ecjLE[NL] , (6.2) 

®The formula can be generalized in a straightforward manner if the efficiencies are different. 
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Figure 5. Expected precision Pz for the SU{5) test on flavour-changing stop decays Eq. (6.3). 
Ny and Nl are the numbers of observed stop decays to bino and wino, respectively. The charm 
fraction is fixed to 7 = 0.5 and the charm-tagging efficiency to Ec = 0.5. Plain (dotted) lines 
denote 2a (Scr) significance respectively. Blue, purple, range and red lines respectively show P 2.3 = 
(100%, 50%, 20%, 10%) isolines of the expected precision. 


where 'yY,L are the actual fractions of charm events. As shown in Ref. [33], the SU{5) 
hypothesis implies = IL ^ 7 - The SU{5) test is therefore 


Ni ~ 


(6.3) 


Let us remark that the large QCD error on the underlying cross-sections roughly cancels 
out in the above ratios of event rates. Moreover, no information on the stop mixing angle 
nor the stop masses is necessary to perform the test. 

Following Sec. 4.4, we define the normalized test quantity as 


R = 


Nl Ni 


/ \Ni Ni) 


such that 


E[R] = 


l7r -7 l| 


lY + 7L 

The expected precision (see Sec. 4.4) associated with this test is found to be 


(6.4) 


(6.5) 


Pz 


2,IP\Ny^NlJ V7 7-1 (7-l)V 


( 6 . 6 ) 


Note that Pz —00 in the limits 7 —)■ 0 (no charm jets expected), 7 —)■ 1 (no up jets 
expected), or ec —?■ 0 (no charm-tagging). The expected precision is shown in Fig. 5, where 
the actual charm fraction is set to 7 = 0.5, and the charm-tagging efficiency is set to 
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ec = 0.5. We see that Ny = > 27 events are sufficient to start to probe the relation at 

3a significance, i.e. to have P 3 < 100%. 

Testing the relation with 50%, 20%, or 10% precision at 3a requires respectively Ny ~ 
Nl ~ 110, 675, or 2700 events. For comparison, assuming stop flavour-violating branching 
ratios of 0.05 and 300 fb“^ of integrated luminosity, we expect about 1340, 130, and 11 
events for stop masses of = 700, 1000, and 1400 GeV, respectively. 


6.2 The 2 ^ case 

For this mass ordering, the stops decay only into the bino according to —)• 

Performing top polarimetry on a decaying stop potentially gives access to the stop mixing 
angle [83]. Here we point out that the same kind of procedure also provides a 5f7(5)-test. 

In the language of the up-squark effective Lagrangian of Sec. 4, the operator coupling 
to is O oc B{tL, — 4:t_R) [37]. The factor —4 comes from the different hypercharges 

of the left and right top components. In the stop effective Lagrangian, the effective mass 
term mixes ii and tR. The mass eigenstates ia, 4 are obtained by 


(I) - (-:::) (t) ^ 

where se and ce are sine and cosine of the stop mixing angle 9t. No implicit assumption is 
made about the mass ordering of ia, it,, i.e. we do not know in principle whether 4 = 4, 
4 = 4 or vice versa The stop angle can be large and is important for phenomenology. 

Using the effective Lagrangian of the stop, it can be shown [33] that the top-stop-bino 
coupling is distorted with a structure 


B{tL,-4tR)JCR{9t) 



K = 



( 6 . 8 ) 


if the SU{5) hypothesis is true. The parameter x encloses the effect of effective operators 
and therefore x <C 1. The matrix JC is instead non-symmetric if the SU{5) hypothesis is 
not satisfied. At leading order in the EFT expansion, i.e. neglecting effective operators, we 
have /C = I 2 . 

As in Sec. 5, let us assume that the spin of the tops is analyzed through distributions 
of the form (1 -|- KPtz) with 2 : G [—1,1]. The decays of the stops 4 and 4 , leading to the 
event rates Na and Nh, are then splitted over the domains = [— 1 , 0 ] and = [ 0 , 1 ], 
such that Na = Na+ + Na- and Nb = Nb+ + Nb-. The expected event rates depend on x 
and 9t, and are given in App. B. 

The stop angle is obtained up to 0{x) corrections using the event rate expressions with 
X = 0 (see App. B), i.e. at leading order in the EFT expansion. Note that only one stop 
is necessary for that purpose. The stop mixing angle is for example obtained from the 4 
forward-backward asymmetry, that satisfies 


nAa] 


K 15 — 17c26) 

2 17 - 15c20 ■ 


®The stop rotation matrix is real as we consider real SUSY breaking terms. 


(6.9) 
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Figure 6. Information content of the event rates Na^b± with respect to the x parameter (see 
Eq. (6.8)). Plain, dashed, dotted and dot-dashed lines correspond to I[Na+], I[Na-], I[Ni,+], and 
I[Ni,-], respectively. The spin-analyzing power is set to k = 1. 


The X parameter is obtained using the full event rates expressions of App. B, i.e. at next-to- 
leading order in the EFT expansion. The measurement of Na^b± provides two inequivalent 
ways of knowing x. The event rates A^a,b± satisfy therefore one non-trivial relation if the 
SU{5) hypothesis is verified. This SU{5) relation can be put in the normalized form 


{Ag - Ab)(1606 + 450c4e) + (1028At - 2040(A^ + Ab))c2e) 

(765 + 255c4e)K + 2040{Ag - Ab)c2e - (1606 + 450c4e)(Aa + Af,) ’ 
where we have introduced the usual asymmetries 


A 


a 


|iVa+ - iVa- 

Na+ + Ng- 


and Ab 


|iVfe+ - Nb- 

Nb+ + Nb- 


( 6 . 10 ) 


( 6 . 11 ) 


One has E[i?] = 0 if the SU{5) hypothesis is satisfied. 

The information about x carried by the expected event rates is conveniently measured 
by I[Ng^b±] = \dlogNi^ 2 ±/dx\. It turns out that this information depends crucially on the 
stop mixing angle. This can be seen in Fig. 6, where we show I[Ng^b±] as a function of the 
latter. The information becomes small for 0* ~ 0 (no stop mixing) and vanishes exactly for 
9t = 7r/4 (maximal stop mixing). In between these two limit cases one has /[iVaj-] = 0(1), 
I[Nb±] = 0(0.1), and the reverse for the interval [vr/d,7r/2]. 

The associated expected precision is found to be approximatively 

p (17-hl5c20) p212 - 739^2 - 1020(^2-4)c2e +(900- 289K2)c40y/2 

^ 255\/2k(3-7c40) V ) 

( 6 . 12 ) 

in the 9t G [0,7r/4] case. We have —)■ oo for k —)■ 0 as the test relies on top polarimetry. 
For 9t G [0,7r/4], the event rates Nb± are less sensitive to x than the event rates Ng±. The 
expected precision depends thus mainly on the amount of tb produced. For this reason we 
drop the l/Ng term is the equation above. The exact formula is given in App. B. 

As for 9t G [0,7r/4] more ib than ig are needed, scenarios where the 4 is the lightest are 
more interesting. For these values of the mixing angle, this lightest stop is mainly right- 
handed. The same reasoning can be applied for 9t G [vr/d, 7r/2]. This time, a larger amount 
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Figure 7. Expected precision Pz for the SU{5) test on polarization of stop decays Eq. (6.10) as a 
function of the spin-analyzer efficiency k and of the amount of observed tb decays Nb- The number 
of ta decays is fixed to Na = 20 (thick lines) and Na = 10000 (thin lines). The stop mixing angle is 
fixed to 6 t = 0.4. Plain (dotted) lines denote 2a {3a) significance respectively. Blue, purple, range 
and red lines respectively show P 2.3 = (100%, 50%, 20%, 10%) isolines of expected precision. 


of ta is necessary. A spectrum where ia is the lightest stop is therefore more favorable for 
the SU{5) test. Again, this lightest stop would be mainly right-handed. We conclude that 
scenarios with a light mainly right-handed stop are always more appropriate to carry out 
this SU{5) test, for any value of the stop angle. 

The expected precision is shown in Fig. 7 for 6 t = 0.4, as well as for Na = 20 and 1000. 
We see that with a spin-analyzer of efficiency k = 0.5 and Na = 20, Nb > 137 events are 
needed to probe the relation at 3cr significance, i.e. to have P 3 < 1. Testing the relation 
with 50% or 20% precision at 3a requires Nb ~ 589 or 7560. For comparison, for 300 fb~^ 
of integrated luminosity, one expects about 26700, 2580, 213, and 24 events for respectively 
rrif = 700, 1000, 1400, and 1800 GeV. 


7 Top-charm supersymmetry 


Let us now focus on classes of models that feature a heavy first generation of up-type 
squarks, i.e. with only the stop and scharm states potentially accessible at the LHC. In the 
super-CKM basis, this means assuming up-squark mass terms of the form 


2 2 

tuq = mu 


'^A 2 0 o\ 
0 0 0 
^ 0 0 oy 


+ 0(M|usy) 


/o 0 0 

0 A22 A32 

yo A23 A33^ 


(7.1) 


where the Xij form in general a hermitian matrix of 0(1) parameters and A ^ Msusy 
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Such a framework is immune against mixing induced from the up-squark 

sector [46] because the first up-squark generation is heavy. We assume that the down-type 
squarks are either aligned or heavy in order to avoid FCNC induced by the down sector. 
Let us recall that, in the 5?7(5)-like GUT context, slepton mass matrices are related to the 
squark mass matrices. 

Phenomenologically, this top-charm SUSY framework constitutes an ideal playground 
to gain knowledge about how to test the SU{5)-like GUT hypothesis of our interest. It is 
also useful in order to carry out stop searches at the LHC taking into account flavour vio¬ 
lation [84, 85]. In particular, large top-charm mixing is found to both improve naturalness 
and to relax the constraints on stop masses [ 86 ]. 

7.1 Extradimensional realizations 

Here we provide some elements of model-building to demonstrate how the top-charm SUSY 
structure can appear in presence of a compact extra-dimension. We focus on the features 
that make the mass matrix Eq. (7.1) appear. The discussion below is not meant to be 
exhaustive, nor fully detailed. Although not all details are discussed, care has been taken 
that the elements proposed are compatible with the observed SM flavour structure and 
that o^-terms are potentially generated. The reader with no interest in model-building can 
safely skip this Subsection. 


• A flat O(TeV) extra-dimension 

Consider a flat extra-dimension y G [ 0 , 7 ri?], with 1/i? of order of few TeVs. Let the 
Higgs Hu and the first generation lOi come from hypermultiplets propagating in the 
bulk, while the two other generations are brane-localized. Supersymmetry is broken 
by twisted boundary conditions (BCs), i.e. Scherk-Schwarz breaking with maximal 
twist (see Refs. [87-104] and the review [105]). The first generation up-squarks get a 
tree-level mass of 1/2R and come within a SUSY N = 2 hypermultiplet. The masses 
of the second and third generations are generated at one-loop [ 88 ], 


mg 


mjj 


diag 


J_ 7C(3) 

4R2’ 167r4R2^’ 


7C(3) ' 


(7.2) 


where C is a model-dependent factor containing group theory invariants. The hier¬ 
archy of Eq. (7.1) is generated due to this loop suppression. 


• A flat O(Mgut) extra-dimension with gauge-Higgs unification 

Let us consider now 1/72 ~ Mqut- The MSSM Higgses come as zero modes of the 
scalar component of the M = 2 vector multiplet. Matter fields are in the bulk with 
some exponential profile controlled by their bulk mass Oj/ttR. The Yukawa couplings 
come from the overlap with matter in the bulk (see Refs. [106-108]). Supersymmetry 
is broken by radion mediation parameterized by the radion F-term Ft (i.e. a Scherk- 
Schwarz with small twist) and by a spurion Fx on the y = 0 brane. The lOi is 
localized towards the y = 0 brane using a oi > 0 bulk mass, IO 2 is localized towards 
the y = ttR brane using 02 < 0, and the IO 3 has flat profile, 03 = 0. A correct Yukawa 
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hierarchy is obtained in that way. For the up-squark masses one obtains 


10 




ml 


diag 


/ Fx V ai 


0 



(Fx 1 
VM2 7ri?M* 


2l 


(7.3) 


where the second-generation term is exponentially suppressed with respect to {Ft/ 2 R)‘^ 
The five-dimensional cutoff M* is such that vriJM* is of order of few units. It is nat¬ 
ural to assume Fx/M^ ~ Ft/2R. In order to reproduce the up-squark Yukawa, one 
has ai ~ 10 , so that there is one order of magnitude between the first and third 
generation. The gaugino having a flat profile, its mass is of same order as the third 
generation of the soft masses ~ Below the GUT scale, the theory re¬ 

duces to the MSSM. The RG flow of the up-squark soft masses is dominated by the 
gluino, whose contribution is universal. The low-energy matrix is roughly given by 
C /3 ~ '^’q u 3 F with Ks = (4.5 —6.5) [37]. One has therefore a low-energy 
matrix of the form Eq. (7.1). 


• A O(Mgut) weakly warped extra-dimension 

Assume a warped extra-dimension with AdS curvature k, with first KK excitations 
iT^KK ~ '^ke at the GUT scale, and the warp factor e ~ 1/20. The MSSM Hig¬ 
gses (possibly containing pNGBs) are localized on the IR brane. The bulk masses 
naturally generate the SM flavour hierarchy by taking 0(1) values [59, 109]. As¬ 
sume that an IR spurion Fx breaks the SM flavour symmetry, giving a soft mass 
0{{Fx/M^)'^) only to the first generation of squarks and to nothing else. Radius 
stabilization and a vanishing cosmological constant [59, 110, 111] typically implies 
that radion mediation is suppressed by a warped factor with respect to the IR brane 
breaking Ft/2R ~ eFx/M^. We therefore have 


2 2 

mq = mu 



2 

diag[l,e^e2]. 


(7.4) 


The gaugino mass is generated by the radion, ~ Ft/2R. The low-energy mass 
matrix has therefore the structure given in Eq. (7.1). 


• A O(TeV) partly-supersymmetric warped extra-dimension 

Finally, let us consider a warped extra-dimension with first KK excitations mxK ~ 
Trke about 100 TeV and the warp factor about e = 10“^®. The Higgs lies on the IR 
brane. Supersymmetry is broken by boundary conditions. Fermion component of 
the hypermultiplet are given (±, ±) BCs, while scalar components have twisted BGs 
(±,^) or (^,±). This setup has a natural holographic interpration has a J\f = 2 
strongly interacting GFT coupled to a non-SUSY elementary sector [80, 112, 113]. 
Bulk masses of 10* are set in order to reproduce the Yukawa hierarchy, with ci > C 2 > 
C 3 ~ 1/2. Assuming left-handed zero modes for the quarks, one sets (^, ±) BGs to the 
1st generation and (±, =f) to the two other generations. From the latter BGs a light 
mode emerges, with m « (4c^ — \')^l‘^ke^^^l’^ for c > 1/2 and m k, V^/j logej^/^/ce for 

do not write the cross-term FtFx for simplicity. 
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c = 1/2. This mode does not exist for the first generation, whose first mode is given 
by mxK = (ci/2 + l)7rA:e. The bulk masses being set in order to reproduce the up 
quark masses pattern, one ends up with 




= ml 


(/se)^ diag 


18, 4-10 


-4 


0.05 


(7.5) 


Notice one obtains an inverted hierarchy where the 2nd generation is the lightest. 
Extra source of SUSY breaking like the universal gravity contributions discovered 
in Ref. [81] could be useful to lift the scharm masses and to generate a-terms. The 
phenomenology of this last model is rather intriguing as the light stops and scharms 
are accompanied by their complex conjugate, as a remnant of the broken A/" = 2 
hypermultiplet. We stress that this model would deserve a dedicated study. 


7.2 Effective Lagrangian 


The effective Lagrangian derived in Sec. 4 applies directly to the top-charm SUSY spectrum, 
identifying the first generation up-squarks as the heavy field 4> = {ul, ur), and the second 
and third generation up-squarks as the light field (j) = {cL,tL,cji,tR). The blocks of the 
mass matrix have then the form 



mfg 


m 


15 


m 


'45 



Furthermore, the low-energy SU{5) relation ~ 




^22 ^"23 


m 


33 


2 

2 \ 

^25 

^26^ 

2 

9 

^35 

COI 

2 

2 

^55 

^56 


mlej 

(7.6) 


of our interest translates into 


77124 


m?5 


and 


mle 


mu- 


(7.7) 


Using also the SU{5) relation Mq ~ (valid only for the two first generations), we have 
in addition 777^2 ~ ^45 and 777^3 ~ TTi^g. It is therefore natural to scrutinize the effects of the 
virtual first generation up-squarks on the light top-charm squarks. The OM~‘^M‘^(j) term 
of the effective Lagrangian Eq. (4.3) induces flavour-changing decays of the light top-charm 
squarks into uB and uW. Distinguishing between the initial c and i seems difficult, and 
we therefore do not pursue this direction. Note that a distorsion of the flavour-conserving 
couplings coming from the 0(1 — term in Eq. (4.3) is also present. 

However, unlike for the Natural SUSY case, looking at the higher-dimensional opera¬ 
tors is not the only possibility available, because a 517(5) information from also 

remains at leading order in the low-energy mass matrix M^. This is the relation of the 
top-charm sector m^Q = m\^. We have therefore the possibility of testing the SU{5) hy¬ 
pothesis using only the real up-squarks. From now on we thus focus only on the top-charm 
sector. 

^^One uses Ci « 0.7, C 2 « 0.6, Ci « 0.5. One has in particular | loge|(2c2 — = ydvt- 
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7.3 SU{5) test through Higgs production 

Let us first consider a case where all stop and scharm masses are nearly degenerate. This 
possibility happens in particular in low-energy GUTs, where no large stop mixing is needed 
in order to have the correct Higgs mass, see e.g. Ref. [104]. Nearly-degenerate squarks imply 
that the mass insertion approximation (MIA) is valid for the sector of stops and scharms, 
and ~ = nriq. These states should be produced in an eqnally 

abundant way at the LHC, as their production occurs mainly through gluon fusion. 

The off-diagonal elements of the up-type trilinear matrix are identihed with mass 
insertions, = Vu{au)ij The SU{5) hypothesis then implies 

= (i„“)32. (7.8) 

To experimentally test this relation, one may scrutinize the composition of the stop and 
scharm eigenstates. At hrst view, even in the MIA, such an analysis seems difficult because 
of the presence of additional mass-insertions ^ 23 ^: (^ 22 ^ 33 . To overcome this issue, 

one should note the fundamental difference between 62 ^’^^ and The former relates 
to a truly bilinear term, i.e. the scalar masses, while the latter is induced by a trilinear 
term, the squark-Higgs scalar coupling. This fact is somehow hidden if one lets the Higgs 
be on its VEV. The physical Higgs exclusively couples to the LR components of the squark 
eigenstates. This can be seen from the complete low-energy mass matrix in Eq. (2.5). In 
order to set up a SU{5) test, one may therefore use the Higgs as a probe of the squark 
eigenstates. Detecting a Higgs gives a direct access to the coupling (<5^^)23, i.e. to ( 0 ^) 23 . 

The LHC SUSY processes of interest are thus stop and scharm pair production, followed 
by a flavour-violating decay into a squark and a Higgs-boson in one of the decay chains. 
These processes are depicted in Fig. 8 . We further assume that the squarks decay into the 
bino. These processes can be identified requiring a single top, a hard jet (from a charm 
quark), a Higgs, and large missing transverse energy (Et). Higgs production through 
up-squark flavour-violating decays have been also studied in Refs. [85]. Note that in the 
degenerate case, not all particles can be on-shell in the decay chain producing the Higgs. 

As in previous cases, this test has to rely on a distinction between the chiralities, which 
is possible only for the top quark. Reconstructing the events is necessary in order to select 
the ones where the Higgs comes from cl,r —>■ htn^L and reject the ones from —)• hcR^R. 

There is in principle enough information from the decay chain to carry out this distinction. 
The former of these processes is shown in first row of Fig. 8 . The latter is shown in second 
row. Other processes leading to the same final states are also possible but are suppressed 
by extra mass-insertions (see second row of Fig. 8 ). 

Provided that the cascade decay with cr —)• hiji can be isolated, top polarimetry then 
readily provides a 5f7(5)-test, as BR(cl —)• hiji) oc and BR(ci? —>■ htL) oc 

Denoting the event number from the relevant flavour-changing Higgs decay chain as 
top polarimetry provides a splitting of the events into N^tc = Nhtc,+ + ^htc- (see Subsec. 
5.1). The SU{5) test then takes the form 

jy \^htc,+ ~ ^htc,— \ 

- AvPnw ^ * ’ 
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Figure 8. Cascade decays in case of a nearly degenerate top-charm spectrum. The dots represent 
mass-insertions. First row: cascade decays used for the SU{5) test defined in Eq. (7.9). Second 
row: other processes leading to the same final state. 


This situation is similar to the one of Sec. 5 and will not be further discussed. 

Instead we focus on a somewhat different spectrum, where stop mixing is large while 
the scharms are nearly degenerate. The MIA applies to the scharm sector, but not inside 
the stop sector. Instead, we rotate exactly the stops into their mass eigenbasis. After 
rotating the stops (see Sec. 6 for the definition of the stop angle), the scharm-stop mass 
matrix takes the form 


( m? 

V 


9 


0 

0 

9 

m- 

Cl 

0 


0 \ 

0 

0 




V 



0 


b^^ce - b^^SQ b^^se + b^^ce \ 
b^^ce - b^^se b^fse + b§^ce 
0 0 

0 0 / 


(7.10) 


Following the MIA approach, the first matrix above corresponds to the squark mass eigen¬ 
values, while the second matrix is treated as a mass insertion. By definition, m? > m? . 
Note that the MIA is expected to be valid to a good precision for ~ 2m\ [49, 50]. 

For the rest of this Section we focus on the case mf > m? > mf . Note that this ordering 

12 r “ 

would happen with degenerate stop-scharm soft masses and a hierarchical On with large 
(3,3) element. Again, this scenario can naturally happen in the five-dimensional GUTs 
discussed in Subsec. 7.1. 
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Figure 9. Cascade decays used for the SU{5) test based on Higgs detection defined in Eq. (7.13). 


The physical Higgs couples only to the left-right mixing terms 6 ^^. The vertices hciii, 
hcLt 2 , hcjiii, hcfii 2 are respectively proportional to |523^S6ip, |(^ 32 ^cep, |<532^seP- 

Within the given mass ordering, the SUSY cascade decays are rather different than from 
the degenerate case discussed above. Flavour changing scharm decays going through t 2 
are now suppressed because of > m?. As a consequence, contrary to the degenerate 
spectrum discussed above, top polarimetry is not useful anymore. On the other hand, real 
decays t 2 —^ he, and c —)■ hii are now opened. 

We require again a single top, a hard jet, a Higgs, and large missing Sy from both 
sides of the decay chains. Two event topologies lead to this final state: the Higgs can 
either come from t 2 —)• or from cl,r —?■ hti. These processes are shown in first 

row of Fig. 9. The first of these diagrams is proportional to -|- |532^pSg, while the 

second is proportional to \^ 2 ^\^s^ + \6^<^\‘^(^. These two types of events can be disentangled 
using the topology of the decay chain. We denote the event rates associated with these two 
diagrams N^j and Nht, respectively. For maximal stop mixing {ce = sq), the two quantities 
become equal so that the power of the test is expected to vanish. 

Contrary to the degenerate case, the stops and scharms have different production rates. 
Moreover the theoretical predictions suffer from a large QCD error. One way to proceed 
is to normalize N^j, N^t using observed data. In order to normalize N^t, we ask for the 
measurement of flavour-conserving decay chains of c-pairs into two jets plus large missing 
E'j'. The corresponding event rate is noted Njj. Because of stop mixing, the same process 
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Figure 10. Expected precision Pz for the SU{5) test based on Higgs detection in SUSY cascade 
decays Eq. (7.13). Nhj and Nht are the number of observed cascade decays from t 2 and cl,r pair 
production, respectively. One fixed the stop angle to 0* = 0.4. Plain (dotted) lines denote 2a {3a) 
significance, respectively. Purple, orange and red lines respectively show P 2.3 = (50%, 20%, 10%) 
isolines of expected precision. 


cannot be used to normalize Nhj. Instead we ask for one of the two t 2 to decay into Zti. 
This event rate is noted Nzt- These processes are depicted in the second row of Fig. 9. 

Normalizing Nhj by Nzt and Nht by Njj cancels the cross-sections, leaving only the 
ratio of partial decay widths, 


E 


r(t2 — hCL,R) 

E 

\Nht] 

_ r(cL,ij —)• hti) 


[l^ztl 

T{t2^Zti) ’ 



^{cl,r cB) 


We now consider the quantity 

_ T{cL,R^hh) Tjh^Zii) _ m? - m| m| - 
“ r(t2 ^ hcL,R) T{cl,r ^cB) ^ m| - m? m| - m| ’ 


(7.11) 


(7.12) 


which is obtained when neglecting the SM masses This factor may be evaluated using 
extra information from kinematic analyses, for example using the kinematic edges of the 
ht, hj, Zt invariant masses. 

The normalized SU{5) test then reads 


which satisfies 


R 


1 

C2e 


Nhj Nht 
Njj "^Nzt 




E[i?] = 


14^-4^1 


(7.13) 


(7.14) 


^^Neglecting the SM masses is made only for illustration purpose, 
the analysis. 


this approximation is not needed for 
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The associated expected precision reads 

Z ( \ 1 \ 

Note that Pz —> oo when 9t —?• vr/d, i.e. the power of the test vanishes in the limit of 
maximal stop mixing as expected. We display only the leading statistical uncertainty that 
comes from the small flavour-changing event rates. 

The expected precision is shown in Fig. 10 for the intermediate value Ot = 0.4. Because 
of the mass ordering of this scenario, one expects N^j < N^t because the t 2 is heavier than 
the scharms and thus produced less abundantly. Assuming Nhj <C Nht and 6 t = 0.4, 
testing the relation with 50%, 20%, or 10% precision at 3 (j requires respectively Nhj > 19, 
116, or 464 events. Roughly twice less events are needed if 6 t = 0. For comparison, 
assuming flavour-violating branching ratios of 0.05, and 300 fb“^ of integrated luminosity, 
one expects about 1340, 130, and 11 events for a squark with respectively rriq = 700, 1000, 
and 1400 GeV. 

Other possibilities of normalization of the Nhj, Nht rates are in principle possible - 
using either observed or theoretical event rates. In any case the approach relies on evalu¬ 
ating the appropriate rj parameter such that Eq. (7.14) is satisfied. The expected precision 
Eq. (7.15) is valid as long as the statistical uncertainty from Nhj and Nht dominates. 

Finally, the analysis of the specific mass orderings carried out above can be taken as 
a guideline to build more global SU{5) tests. In particular the flavour-changing squark 
decays with Higgs production are expected to always carry relevant information regarding 
Qu ~ On) for any mass ordering. Assuming an arbitrary mass ordering, a global SU{5) 
test can be setup by putting into a likelihood function the Nhj , Nht event rates plus some 
information to normalize them (i.e. all the potentially useful information found above). 
This global test constitutes a natural extension to the cases analysed above and is subject 
to future work. 

8 Summary and conclusion 

A fairly intriguing feature of the Standard Model (SM) is that quarks and leptons fit 
naturally into complete representations of the SU{5) gauge group. This might be taken as 
a hint that the SM is the low-energy effective theory of a Grand-Unified Theory (GUT), 
either 517(5)-symmetric or symmetric under a group broken to SU{5) at an intermediate 
energy scale. Many realizations of S'C/(5)-like GUTs are possible, and imply a variety 
of low-energy features. Moreover, the large quantum corrections occurring when flowing 
down the SM effective theory from the GUT to the weak scale potentially wipe out a 
lot of precious information. Testing whether the world is microscopically 5C/(5)-like is 
thus a rather challenging task. The large model-dependence coming from both the GUT 
realization and from quantum corrections constitute a irreducible theoretical uncertainty, 
from the point of view of testing the SU{5) hypothesis. 

Low-energy supersymmetry is a natural ingredient of S'[/(5)-like GUTs, as it typically 
favours the high-energy unihcation of gauge couplings. As a happy coincidence, the super¬ 
partners potentially carry information about the possible GUT group. Low-energy SUSY 
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is therefore both favoured by 5C/(5)-hke GUTs and useful to reveal their nature. Still, 
testing the SU{5) hypothesis rises the double challenge of taming the model-dependence 
and finding observations that are doable in the real world. While such a program may seem 
at first view insurmountable, it happens that a new SU{5) relation between the up-squark 
trilinear couplings au = first identified in Ref. [33], opens a number of possibilities. 
This relation is exact at the GUT scale, and is immune to any GUT threshold corrections. 

Moreover, the relation ~ is found to survive through the MSSM renormalization 
group flow, such that it is spoiled by typically 0(1%) of relative error at the TeV scale. 
This relation can be thus taken as a window on the GUT physics, that can be exploited 
to test the 517(5) hypothesis. Because of the peculiar structure of a(j, all the 517(5)- 
tests coming from this relation have to rely on either flavour-violation or chirality flip in 
the sector of up-squarks and quarks. The discussed 517(5)-tests are set up for the Large 
Hadron Gollider (LHG). 

Even though the 517(5) relation is immune to many corrections, a model- 

dependence remains in the up-squark spectrum. Indeed, although certain patterns of the 
up-squark mass matrix can be justified in a number of ways, overall a lot of freedom 
remains. We consider therefore various typical scenarios for the up-squark spectrum. The 
scenarios that we have investigated are “Heavy supersymmetry” (i.e. Msusy TeV, Sec. 
5), “Natural supersymmetry” (i.e. real stops are accessible at the LHC, Sec. 6), and “Top- 
charm supersymmetry” (i.e. both scharm and stops are accessible at the LHC, Sec. 7). 

Depending on the up-squark spectrum, the 517(5) tests involve only virtual squarks, 
only real squarks, or both real and virtual squarks. Whenever the virtual squarks are heav¬ 
ier than the TeV scale, they can be integrated out to form a TeV-scale effective Lagrangian. 
The effect of virtual squarks is then enclosed into effective operators, either generated at 
tree-level or at one-loop. A complementary expansion is the mass insertion approximation, 
that we use both in the broken and unbroken EW phase. 

In many cases, the 517(5) tests that we find consist in determining whether a relation 
among certain observables is satisfied or not. In order to quantify the feasibility of a test, 
we introduce a systematic procedure relying on the frequentist p-value. The associated 
expected precision tells, for a given amount of data, up to which magnitude a violation of 
the 517(5) relation can be assessed within a given statistical significance. 

In the “Heavy SUSY” case, it appears that the 517(5) hypothesis can be tested using 
flavour-changing dipole operators in the up-sector. These dipoles from SUSY are loop¬ 
generated. A top quark has to be present in the dipole in order to use top-polarimetry. 
The expected sensitivity as a function of signal and background event rates is evaluated. 
The common LHC searches from either anomalous single top production or anomalous 
flavour changing top decays seem however not sensitive enough to be able to discover 
SUSY dipoles. 

We propose thus an alternative kind of precision measurements relying on ultraperiph¬ 
eral collisions between proton and/or heavy ions. A set of exclusive processes is proposed 
in order to probe the SUSY dipoles operators. We also point out that, contrary to a seem¬ 
ingly widespread belief, the W / Z-boson fluxes from proton and heavy ions are coherent in 
the UPC regime. The vector boson fusion cross-sections can thus be expected to be rather 
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large. 

In the “Natural SUSY” scenario, already investigated in Ref. [33], the two mass order¬ 
ings mi^ 2 ^ ^BW ^ ^ii 2 ^ considered. A test associated to the former 

relies on flavour-changing stop decays and on charm-tagging, while no knowledge of stop 
masses nor the stop angle is needed. We evaluate the expected precision as a function of 
charm-tagging efficiency and event rates. 

For the second mass ordering, the SU{5) test relies on top polarimetry. Its potential 
greatly depends on the stop mixing angle, and vanishes in the limit of no stop mixing 
and maximal stop mixing. In contrast, the potential for arbitrary stop mixing is more 
promising. We evaluate the expected precision as a function of the stop mixing angle, 
event rates and spin-analyzing power. A few events from the mainly left-handed stop are 
sufficient, while a more copious production is needed for the mainly right-handed stop. 
This test is favoured by a SUSY spectrum with a mainly right-handed lightest stop. 

We finally introduce the “Top-charm SUSY” scenario, featuring a heavy first genera¬ 
tion of up-squarks, and demonstrate that it can arise in a variety of five-dimensional GUT 
realizations. SU (5) tests can be found using the SM Higgs boson to probe the LR compo¬ 
nents of the squarks. The tests therefore rely on cascade decays involving Higgs production 
and flavour-violation. For degenerate scharms and stops, a test using top polarimetry is 
presented. For the case of > ruci r > ^ test relying on the reconstruction of the 

decay chains involving the Higgs is discussed. The expected precision as a function of the 
stop angle and event rates is evaluated. 

The SU{5) tests that appear in the various SUSY scenarios considered are summarized 
in Tab. 4. The typical amount of events needed to reach an expected precision of 50% at 
3a is also shown for each of the tests. For these numbers the experimental conditions are 
chosen to be the most favorable, e.g k = 1. The number of needed event ranges roughly 
from about 10 to about 100. For the required precision, the tests involving top polarimetry 
require typically 0(100) events, because they rely on the shape of kinematic distributions. 
All these tests rely on the existence of a given relation among certain observables. In cases 
where no such relation is available, a more global hypothesis testing has to be done. This 
will be the subject of a further work [114]. 
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Heavy 

Natural SUSY 

Top-charm SUSY 


SUSY 


> bUb 

mr ~ rricj „ 

B^i2>^CL.R 

> "iq 

Squarks involved 

virtual 

virtual/real 

real 

Top polarimetry 

yes 

no 

yes 

yes 

no 

Charm-tagging 

no 

yes 

no 

no 

no 

Higgs detection 

no 

no 

no 

yes 

yes 

0t-dependence 

no 

no 

yes 

no 

yes 

Fg = 50% 

144 

72 

108 

144 

10 


Table 4. Summary of the S'C/(5)-tests appearing in the various SUSY scenarios considered. The 
last line shows the typical number of events needed to reach a 50% precision at Scr. 
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dipole form factors 

factors from penguin diagrams are [57, 58] 
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21(1 -x)4 7(1 -x)5 

12(ll + x) 6(9 + 16 x-x2) 


5(1-x)3 5(1-x)4 

10 (1 - 8x - 17x2) _ 20 x 2(3 + x) 
3(1 — x)^ (1 — x)® 

6(l + 5x) 12x(2 + x) 


logx, 

logx, 


(1 — x)" 


+ 


(1 — xY 


■ log X. 
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B Event rates for ta^ tB 

The expected event rates have the form 


.r Ar l-«/2 1 +k/2 
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(B.l) 
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(B.3) 

I + k/2 I-k/2 

Nb- = NbL 2 ^’^B 2 ’ 
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NbL = Nh^, 
Nm = Nb 


16 (s 0 - xcef 
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Here Na^h is the production rate of tafi- The expected precision associated with the test of 
Eq. (6.10) reads 


Pz = 


255\/2 K (3 + C4e) . 


(17 - 15c2e)2(3212 - 739^2 + 1020(«;2 - 4)c20 + (900 - 289K^)c^g) 
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(17 + 15c2e)2(3212 - 739^2 _ 1020(^2 - 4)c20 + (900 - 289R?)cie)^ 
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